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Abstract 

We study second order phase transitions in non-conformal holographic models of 
gauge theory/string theory correspondence at finite temperature and zero chemical 
potential. We compute critical exponents of the bulk viscosity near the transition 
and interpret our results in the framework of available models of dynamical critical 
phenomena. Intriguingly, although some of the models we discuss belong to different 
static universality classes, they appear to share the same dynamical critical exponent. 
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1 Introduction 



Gauge theory/string theory correspondence of Maldacena [1] presents a useful theo- 
retical framework to study strong coupling dynamics of gauge theories. In this paper 
we use it to study bulk viscosity of non-conformal gauge theory plasma at criticality 
of a second order phase transition. We focus on finite temperature phase transitions 
at zero chemical potential^. Furthermore, since the models we consider do not have 
long-range (Goldstone boson) modes at the criticality, their hydrodynamics is that of 
a standard viscous relativistic fluid , i.e, the local stress-energy tensor is given by 

T^" ^eu^u" + P(e) A'^'^ - r]{e)a^'' - C(e) A'^" (V • u) , 

2 (1-1) 

where e and P{e) are the local energy density and pressure, is the local d-velocity of 
the plasma, and 77(e) and C(e) are the shear and the bulk viscosities correspondingly. A 
plasma with the stress-energy tensor (1.1) allows for a propagation of the hydrodynamic 
sound waves with the following dispersion relation 

It) = c2q-zrq2 + 0(q3) , (1.2) 
where Cg is the speed of sound and F is the sound wave attenuation, 

r = 2.5(14 + ^), (1.3) 



de J J, s \d — 1 2r] 

and Vo = uj/{2ttT) and q = |g|/(27rT). Thus, if the dispersion relation for the sound 
waves (1.2) in plasma can be determined from first principles, the bulk viscosity can be 
computed from (1.3)^. A holographic correspondence of Maldacena provides a frame- 
work for such first principle computations [3]: the dispersion relation of the sound 
waves in strongly coupled plasma is identified with the dispersion relation of a hy- 
drodynamic graviton quasinormal mode (of a certain polarization) in gravitational 
background holographically dual to the equilibrium thermal state of the plasma. 
In a conformal plasma the trace of the stress-energy tensor vanishes. As a result, 

CFT "J- '/ CFT 



"'^ Critical phenomena in strongly coupled plasma at finite density will be discussed in [2]. 
^One also needs a ratio of shear viscosity to the entropy density. In models we analyze this ratio 
is universal [5, 6, 7]. 
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However, once the scale invariance is explicitly broken, the bulk viscosity is generically 
nonzero. While for most physical substances the bulk viscosity is smaller or at most 
of the same order as the shear viscosity, ^He in the vicinity of the critical liquid- vapor 
point exhibits the ratio of the bulk-to-shear viscosities in the excess of a million [4]. 
Below, we present the first example of a holographic model with divergent ratio of ^ 
in the vicinity of a second order phase transition. 

Techniques for the holographic computations of the bulk viscosity in non-conformal 
plasma were developed in [8]. In what follows we omit all the technical details of the 
analysis. In the next section we mention three different proposals for the scaling of 
bulk viscosity at the criticality of a second order phase transition [9, 10, 11]. In section 
3, building up on the previous work [12, 13, 14], we identify the second order phase 
transitions (and the static universality classes) in holographic models of gauge-gravity 
correspondence: the J\f = 2* gauge theory plasma [15, 16, 17], the finite temperature 
Klebanov-Tscytlin (cascading) gauge theory plasma [18, 19, 20, 21, 22, 23], and the 
phase transition in "exotic black holes" wc proposed earlier [24]. In section 4 we present 
results for the scaling of bulk viscosity at the criticality in these models. While the 
computations of bulk viscosity in J\f — 2* and the cascading gauge theories were done 
previously (see [8, 13, 25, 26, 14] ), our result for the bulk viscosity of plasma dual to 
"exotic black holes" of [24] is new. We conclude in section 5 by interpreting the scaling 
of the bulk viscosity in holographic models in the framework of models of dynamical 
critical phenomena [9, 10, 11]. 

2 Models of bulk viscosity at criticality 

We end up mapping the second order phase transitions in holographic models of 
gauge/gravity correspondence to the ferromagnetic phase transition inp = d — 1 spatial 
dimensions — thus we use the nomenclature of the latter. 

Consider the Gibbs free energy density W = W{T, Ti) which is the difference be- 
tween the free energy densities of the ordered and the disordered phase 

>V(T, H) = n,{T, H) - QdiT, H) (2.1) 

as a function of temperature T and the (generahzed) external magnetic field Ti.. The 
(generalized) spontaneous magnetization Ai determines the response of the free energy 
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to the changes in the external control parameter, as 

--mi- 

The energy density e, the entropy density s, and the spontaneous magnetization M. 
satisfy the basic thermodynamic relation 

W = e-sT -MH, (2.3) 

with the first law of thermodynamics 

dW^-sdT-Mdn. (2.4) 

For convenience, we introduce a reduced temperature t 

(2.5) 

where T^. is the critical temperature of the second order phase transition. For a standard 
ferromagnetic phase transition^ 

t<0: w{t,n)<o, M{t,n = o)^o, ^ 
t>0: w{t,n)>o, Mit,n = o) = Q. 

At a second order phase transition the first derivatives of W are continuous while the 
higher derivatives are not. Under the static scaling hypothesis we have 

w(t, n) = x-p w (xy^t, x^^n) , (2.7) 

for the free energy, and 

G{q, t, n) = A^^^-f G{Xq, xy^t, x^^n) , (2.8) 

for the Fourier transform of the equilibrium two-point correlation function of the mag- 
netization 

G«HA<«A.(0)>oc^^^|^. (2.9) 

The static critical exponents {a, /?, 7, 5, i/, r/} are defined as 

1/3 



cw = -r(^l <x\t\-\ Mcx\t\^, 



(2.10) 



^Some of the holographic phase transitions we discuss are "not standard", see [24]. 
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and 

( g-|r1/C ^ / Q 

G{r)(x{ ' , with (2.11) 

[ \f\-p+^-\ i = 

where c-^ is the specific heat, xr is the magnetic susceptibihty, and ^ is the correla- 
tion length. Given (2.4) and the scaling hypothesis (2.7) and (2.8), the two critical 
exponents {yr, yn} determine the critical exponents (2.10) which identify the standard 
static universality classes: 



a = 2 , /3= , 7 = 

Vt yr yr 

5 = -J2L_ ^ = — , r]=p-2yn + 2. 

p-yn Vt 



(2.12) 



Note that (2.12) implies the following scahng relations 

« + 2/3 + 7 = 2, 7 = /3((^ - 1) = i^(2 - 7?) , 2-a^iyp. (2.13) 

The scaling relations (2.13) can be violated whenever the single-scale hypotheses (2.7) 
and (2.8) break down. In what follows we assume that this is not the case in our 
models. 

A finer classification of the universality classes of systems undergoing second order 
phase transition occurs once the dynamical properties of the system (such as hydro- 
dynamics) are under consideration [27]. Here, the main assumption is that a system 
perturbed away from the equilibrium will relax with a characteristic time scale Tg, 
which diverges near the phase transition with a dynamical critical exponent z: 

Tg oc r oc |^|~'^^ (2.14) 

Correspondingly, the near-equilibrium scaling of the magnetization correlation function 
(2.8) is modified according to 

g{lj, q, t, n) = x^y^--^^' G{X'tj, Xq, xy^t, x^^n) . (2.15) 

In this paper we will be interested in the hydrodynamic aspects of the critical phe- 
nomena, specifically in the scaling of the bulk viscosity at the transition. Although 
second order transitions imply spatial scale invariance (2.11), the latter does not pre- 
clude a non-zero bulk viscosity which necessitates the space-time scale invariance, see 
(1.4). A large bulk viscosity at the transition signals strong coupling between the 
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dilatational modes of the system and its internal degrees of freedom. This is a very 
interesting phenomena given that hydrodynamics is often regarded as an effective low- 
energy description of a system. Without going into details we now mention three 
proposals (Models A,B,C) for the scaling of the bulk viscosity at criticality^. 
■ (Model A): In [9] the authors proposed that the critical behavior of the bulk viscosity 
is governed by the critical exponent of the specific heat 

C oc oc (2.16) 



■ (Model B): In the quasi-particle model under the relaxation time approximation the 
bulk viscosity was argued to scale as [10] 



C oc (2.17) 



■ (Model C): In [11] Onuki argued that the bulk viscosity near the gas-liquid critical 
point (model H in classification of [27]) scales as 

Coc|ir^"+". (2.18) 

Below, we discuss second order phase transitions in holographic models of gauge/gravity 
correspondence. We show that only the scaling (2.18) passes these holographic tests. 

3 Holographic second order phase transitions 

A second order phase transition in a holographic model of gauge/gravity correspon- 
dence was first reported in [28]. This phase transition is associated with the sponta- 
neous breaking of U{1) symmetry, and as a result the hydrodynamics of the model 
must include a long-range Goldstone boson mode. Another type of a holographic sec- 
ond order phase transition was discussed in [29]. The phase transitions in [28] and [29] 
occur both at finite temperature and finite chemical potential. 

We discuss bulk viscosity in second order phase transitions at finite density in [2]. 
Here, we focus on phase transitions at vanishing chemical potential for the conserved 
charges. An example of such transition was reported in [24] (we review the main 

^Hcrc wc indicate only the scaling of the potentially singular part of the bulk viscosity. A finite 
contribution is always implicit. 
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Figure 1: (Colour online) The speed of sound (left plot) and the reduced temperature 
^ (right plot) of the strongly coupled Af = 2* plasma with mj = and 7^ as 
a function of the dual gravitation parameter pn. We identify the blue curve with the 
"ordered" phase and the red curve with the "disordered" phase, see (2.1). 

features of this transition in section 3.3). We show that both the M = 2* gauge theory 
and the cascading gauge theory plasma undergo the second order phase transitions 
which are in the same static universality class. The latter universality class includes 
the model of [29], but is different from the universality class of the model [24]. The 
static universality classes of the holographic models we discuss below are not of the 
mean-field type since the critical exponent 77 (see (2.11)) is nonzero. 

3.1 = 2* plasma 

The holographic renormalization and the thermodynamics of the M = 2* gauge theory 
plasma in the planar limit and at (infinitely) large 't Hooft coupling was discussed pre- 
viously [30, 31, 12]^. This gauge theory is obtained by giving a mass to one of the three 
M = 2 hypermultiplets of A/" = 4 SU{N) supersymmetric Yang-Mills theory. At finite 
temperature the supersymmetry is broken, and one can study the theory with differ- 
ent masses for bosonic and fermionic components of the massive hypermultiplet. In 
[12] the different mass deformations were called the "bosonic" m;,, and the "fermionic" 
TUf. For the range of parameters {T,mh,mf) studied in [12, 13] the Af = 2* plasma is 
always in deconfined phase^. It was further shown that whenever m'j < ml the theory 

'^Thermodynamics of weakly coupled JV = 2 SU (2) Yang-Mills theory was recently discussed in 
[32]. 

^An interesting phase transition in this model was conjectured in [30], see also [32]. 
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Figure 2: (Colour online) Free energy densities flo of the "ordered" phase (blue curves) 
and Qd of the "disordered" phase (red curves) as a function of pu (left plot) and ^ 
(right plot) of the M = 2* plasma with = 0. 



undergoes a phase transition with the vanishing speed of sound. Although the critical 

m? 

temperature T,, of the transition depends on the ratio the general thermodynamic 
and the hydrodynamic features of the model are the same [25]. Thus, without the 
loss of generality we consider the mass deformation rrif = and rrif, ^ 0. Here, the 
transition occurs at [12] 

^^ 2.32591. (3.1) 

We now identify this phase transition as a second order phase transition with the static 
critical exponents 

{a,(3,^,6,u,v)= (^^,^,^,2,^,1^ . (3.2) 

It is convenient to present the thermodynamic data as a function of the dual grav- 
itational parameter pu. For T ^ nib we have [30] 

for all other temperatures see section 3.3 of [12]. The left plot on Fig. 1 represents 
the square of the speed of sound cj. as a function of pn while the right plot on Fig. 1 
represents the reduced temperature ^ as a function of pu. The critical — the lowest — 
temperature corresponds to the maximum on the plot. The plots on Fig. 2 represent the 
(reduced) free energy of the system as a function of pu (left plot) and the reduced 
temperature (right plot). For a given temperature T > Tc there are two phases of 
M = 2* plasma denoted by a blue/red dot. From Fig. 2 it is clear that the transition 
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between the two phases is a continuous one. The "red" phase is perturbatively unstable 



as it has c? 



< — it is natural to identify it with the disordered phase of the effective 



red 



ferromagnet in section 2. Perturbatively stable "blue" phase is thus an ordered one. 
Prom the right plot in Pig. 2 we see that the ordered phase is thermodynamically 
favorable since 

W = - = Q^'"^ - fi^^'^ < , (3.4) 

Note that in this system, contrary to the standard ferromagnet (2.6), the dimensionless 
temperature, defined by (2.5), t > always. Let's denote by the critical value of 
pii, i.e., 

(3.5) 
(3.6) 



T 



pii=Pii 



rUb 



Introducing 

from Pig. 1 it is clear that 
t oc (Apii)', 

Thus, 



oc 



blue 



OC |Apii| oc 



red 



blue 



oc C 



-2 



red 



blue 



OC t' 



-1/2 



(3.7) 



(3. 



red 



where we used the fact that the entropy density s is continuous across the phase 
transition. Comparing (3.8) and (2.10) we determine the critical exponent a [13] 



"=2- 



(3.9) 



To determine the critical exponent (3 we need to identify the control parameter corre- 
sponding to the external magnetic field H of the effective ferromagnet of section 2. We 
propose to identify 

n^rub. (3.10) 

Given (3.1), it follows from (3.7) that 

1 



dn oc -dt oc 



Ap 



11 



Apii 



(3.11) 



The left plot on Pig. 2 implies that 



(3.12) 
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0.5376 




Figure 3: (Colour online) The reduced temperature ^ (left plot) and the free energy 
densities Qo of the "ordered" phase (blue curve, right plot) and Qd of the "disordered" 
phase (red curve, right plot) , see (2.1), of the strongly coupled cascading plasma as a 
function of the dual gravitational parameter kg. 



at the very least. Actually, the best fit to fl shows that^ 



W oc -|Apii 



Thus, using (3.11) and (3.13) 
M 



1 



oc 



dn J Apn 
Comparing (3.14) and (2.10) we conclude that 



/3 = -. 
^ 2 



^ApiiW OC -|Apii| oc -t 



1/2 



(3.13) 



(3.14) 



(3.15) 



The rest of the critical exponents in (3.2) follow from the scaling relations (2.13). 

The second order phase transition in A/" = 2* plasma we just exhibit is remarkably 
similar to the second order phase transition in A/" = 4 supersymmetric Yang-Mills 
plasma with a single U{1) C SU{A)fi R-symmetry chemical potential discussed in [29]. 
In [29] this chemical potential was identified with the external magnetic field 7i of the 
effective ferromagnet of section 2. 



3.2 Cascading gauge theory plasma 

The holographic renormalization and the thermodynamics of the cascading gauge the- 
ory plasma in the planar limit and at (infinitely) large 't Hooft coupling was discussed 

^This is not surprising since fl on the left plot of Fig. 2 can not be an even function of Apu. 
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previously in [20, 21, 22, 33, 23, 14]. The cascading gauge theory [19] is best thought 
of as jV = 1 supersymmetric SU {K + P) x SU {K) gauge theory where the effective 
number of colours K is not constant along the renormalization group flow, but changes 
with energy according to [20, 34, 35] 

K ^ K{n) r^2P^\nf^, (3.16) 

at least when the energy scale /i is much larger than the strong coupling scale A of the 
gauge theory. At zero temperature, and when K{fi) is a multiple of P, cascading gauge 
theory confines in the infrared with the spontaneous breaking of chiral symmetry. The 
high temperature phase of the theory is expected to be that of the deconfined chirally 
symmetric plasma [20]. In was shown in [23] that at 

T = refinement = 0.6141111(3)A (3.17) 

the deconfined chirally symmetric phase undergoes a first-order confinement phase 
transition, with spontaneous breaking of chiral symmetry. For 

Tc 0.8749(0) X T^onfinement ^ T <C T confinement i (3.18) 

the deconfined phase of the cascading plasma, although non-perturbatively unstable 
due to the nucleation of bubbles of the confined phase, remains perturbatively (and 
thermodynamically) stable [14]. In the vicinity of Tc the thermodynamics of the cas- 
cading gauge theory plasma is identical to that of the J\f —2* plasma. Fig. 3 presents 
the reduced temperature ^ (left plot) and the free energy density of the cascading 
plasma as a function of the dual gravitational parameter ks- We use the same colour 
coding as in Fig. 1 and Fig. 2: the "blue" dots indicate stable (ordered) phase and the 
"red" dots indicate unstable (disordered) phase. As in the case ol M — 2* plasma, in 
the ordered (disordered) phase the square of the speed of sound is positive (negative) 
[14]. Once we identify the external magnetic field V, of the effective ferromagnet of 
section 2 as 

7i = A, (3.19) 

we can literally repeat the arguments of section 3.1 and arrive at the same set of static 
critical exponents (3.2). 

3.3 "Exotic black holes" 

Following Gubscr's suggestion [36], we constructed an exotic model of the second or- 
der phase transition in o? = 3 at finite temperature and zero chemical potentials in 
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Figure 4: (Colour online) The free energy densities Qq of the "ordered" phase (red 
curve, left plot) and Qd of the "disordered" phase (purple curve, left plot) as a function 
of the reduced temperature ^ in gauge theory plasma dual to holographic RG flow in 
[24]. The right plot represents the square of (Oi) (which we use an an order parameter 
for the transition) as a function of reduced temperature. The dashed green line is a 
linear fit to (C^) 2. 

[24]. Specifically, we considered relativistic conformal field theory in 2+1 dimensions, 
deformed by a relevant operator Or'- 

HcFT —^H = Hqft + Aj.0r • (3.20) 

Such a deformation softly breaks the scale invariance and induces the renormalization 
group flow. We further assumed that the deformed theory H has an irrelevant op- 
erator Oi that mixes along the RG flow with Or- In the explicit holographic model 
realizing this scenario [24], the irrelevant operator Oi developed a vacuum expectation 
value in the high temperature phase, i.e., for T > T^, spontaneously breaking a dis- 
crete Z2 symmetry of the model. The unusual part of this phase transition was that 
the symmetry broken phase occurs at high temperatures (rather then the low temper- 
atures) and that the broken phase (although perturbatively and thermodynamically 
stable) has higher free energy density than the unbroken phase with (Oi) = 0. The 
various high temperature phases which spontaneously break Z2 symmetry differ in the 
holographic wavefunction of the condensate and their free energies [24]. Without the 
loss of generality we focus here on the model with dim[(9r] = 2, and the broken phase 
of the lowest free energy. Fig. 4 represents the energy densities of the ordered (red 
line, left plot) and the disordered (purple line, left plot) phases as a function of the 
reduced temperature ^. The expectation value of Oi is the order parameter of the 

-L C 
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Figure 5: (Colour online) The speed of sound Cg in gauge theory plasma dual to holo- 
graphic RG flow in [24] as a function of the reduced temperature ^. 

phase transition — the spontaneous magnetization Ai of the effective ferromagnet in 
section 2. An excellent linear fit (the dashed green line) to (Oi)^ (see the right plot on 
Fig. 4) implies that in this model 

M = (Oi) oc |t|^/2^ (3.21) 

resulting in the critical exponent 



P 



(3.22) 



From Fig. 5 it is clear that the speed of sound is finite at the transition^, thus 

cn oc oc (3.23) 

resulting in the critical exponent 

a = 0. (3.24) 
The remaining critical exponents can be obtained from the scaling relations (2.13): 

(3.25) 

Note that the universality classes (3.2) and (3.25) are different. 



(q;,/3,7,5, z/,!]) = ( 0,-, 1,3, 1,1 



4 Holographic bulk viscosity at criticality 

The technique for computing the bulk viscosity in non-conformal holographic models 
from the dispersion relation of the sound waves was developed in [8]^. It was used to 

^Notice that in the disordered phase the speed of sound violates the bound proposed in [37, 38]. 
much simpler computational framework was recently proposed in [39, 40]. 
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compute the bulk viscosity oi N — 2* plasma in [8, 13, 25], and the bulk viscosity of 
the cascading gauge theory plasma in [26, 14]. It is straightforward to generalize the 
method of [8] to the computation of the bulk viscosity of the gauge theory plasma dual 
to the holographic RG flow [24]. The details of the latter analysis will appear elsewhere 
[41], and here we report only the results. 

In the rest of this section for each of the holographic models with a second or- 
der phase transition discussed above we analyze the scaling of the bulk viscosity at 
criticality in the framework of Models A,B,C — (2.16)-(2.18). 

4.1 = 2* plasma 

The main result of [13, 25] was that the ratio of the bulk-to-shear viscosities in jV = 2* 
plasma remains finite across the second order phase transition^'^: 

- oc (4.1) 

Since in holographic phase transitions under the consideration both the entropy density 
and the shear viscosity are finite, we conclude from (4.1) that the bulk viscosity stays 
finite. Thus: 

■ Model A is inconsistent with holographic analysis as it predicts divergent bulk vis- 
cosity, C oc Itl^"*"/^ ; 

■ Model B does not contradict our holographic analysis as it predicts that (.singular oc 

■ Model C agrees with holographic analysis provided the dynamical exponent z is 

z^l. (4.2) 

4.2 Cascading gauge theory plasma 

It was argued in [14] that the bulk viscosity of the cascading gauge theory plasma 
remains finite across the second order phase transition much like in the case of A/" = 2* 
plasma. Thus, cascading gauge theory plasma appears to share not only the static 
universality class with J\f = 2* model, but also the dynamical one — in particular, the 
critical exponent z is given by (4.2). 
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Figure 6: (Colour online) The ratio of bulk-to-sliear viscosities in gauge theory plasma 
dual to the RG flow in [24]. The red curve corresponds to the "ordered" phase of the 
theory, and the purple one to the "disordered" phase. The dashed blue hne indicates 
the bulk viscosity bound proposed in [13]. The dashed green line corresponds to the 
high-temperature approximation to the viscosity ratio given by (4.4). 

4.3 "Exotic black holes" 

Fig. 6 represents the ratio of the bulk-to-shear viscosities of the gauge theory plasma 
dual to the RG flow in [24]. The red curve corresponds to the "ordered" (symmetric) 
phase that the purple curve corresponds to the "disordered" (broken) phase. The 
dashed blue line represents the bulk viscosity bound 



proposed in [13]. The dashed green line represents the bulk-to-shear viscosity ratio for 
a symmetric phase at high temperatures: 



ordered 



V3\ 



2 




cl]+0{{--cl) 1 . (4.4) 



Notice that the bulk viscosity bound is satisfied both in the symmetric and in the 
broken phases. In fact, in the disorder phase (purple curve) the bound is satisfied 
trivially as here > |, see also Fig. 5. 

Fig. 6 shows a rapid rise of the bulk-to-shear viscosity ratio in the disordered (purple 
curve) phase as one approaches the transition. In fact, a detailed analysis presented in 
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The second-order transport coefRcient is nonetheless divergent [42]. 
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Figure 7: (Colour online) The ratio of bulk-to-shear viscosities in gauge theory plasma 
dual to the RG flow in [24] in the vicinity of the second order phase transition as a 
function of the order parameter. The dashed green line represents the linear fit to the 
log-log data plot with the slope of (—1.9999(6)). 



Fig. 7 strongly suggests that bulk viscosity diverges in this case. There, we plot In 
versus In A dashed green line given by 

y = -1.9999(6) a; - 5.7684(5) 
provides an excellent linear fit to the data. Thus, we expect 

oc m 



oc \t\-' 



(4.5) 



(4.6) 



disordered 



where we used (3.21). From (4.6) we conclude: 

■ once again. Model A is inconsistent with our holographic analysis as it predicts a 
finite bulk viscosity at the transition, ( oc ; 

■ Model B is inconsistent as well, as it predicts (singular oc |t|^ ; 

■ Model C agrees with holographic analysis provided the dynamical exponent z is 



1. 



(4.7) 



Rather intriguingly, even though the static universality classes of the M = 2* (and the 
cascading) gauge theory plasma and the exotic plasma discussed in this section are 
different, when the bulk viscosity is interpreted in the framework of Model C, both 
classes have the same dynamical critical exponent z. 
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5 Conclusion 



In this paper we studied bulk viscosity at criticality of the second order phase tran- 
sition in non-conformal models of gauge theory/string theory correspondence at finite 
temperature, but at zero chemical potential for the conserved charges. We identified 
the second order phase transitions in Af — 2* and the cascading gauge theory plasma 
and showed that the two models are in the same static universality class. Somewhat 
surprisingly, this universality class is not of the mean-field theory type. Incidentally, 
the same universality class contains a second order phase transition in A/" = 4 su- 
persymmetric Yang-Mills plasma with a single U{1) C S'f/(4)/j R-symmetry chemical 
potential discussed in [29]. The static universality class of the M = 2* plasma is differ- 
ent from the static universality class of the gauge theory dual to the 3+1 dimensional 
asymptotically ^^154 RG flow discussed in [24]. Having different static universality 
classes in our holographic laboratory allowed us to test different proposals for the scal- 
ing of bulk viscosity in the vicinity of the second order phase transition. We found 
that neither the model of [9] or [10] conforms to the tests. If dynamical scaling of bulk 
viscosity proposed by Onuki [11] is correct, then all our holographic models must have 
the same dynamical critical exponent z = 1. 

In the future, it is important to explicitly verify the scaling relations (2.13) in our 
holographic models. It would be extremely interesting to subject Onuki's theory [11] to 
a direct holographic test by independently computing the dynamical critical exponent 
z from the scaling relation (2.14). Finally, it is interesting to understand the dynamical 
universality class of the holographic models. We already know that in such theories 
the ratio of shear viscosity to the entropy density is universal. Is it possible that the 
universality of the holographic transport at critically is more robust? 
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